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1.  Introduction 


In  this  work  the  full  wave  approach  is  used  to  determine  the 
scattering  cross  sections  for  composite  models  of  non-Gaussian  rough 
surfaces.  In  particular,  it  is  assumed  here  that  the  rough  surface  is 
characterized  by  a  family  of  joint  height  probability  densities  that 
have  been  developed  by  Beckmann  (197 3a, b)  for  non-Gaussian  surfaces. 

These  joint  height  probability  densities  are  expressed  as  an  infinite 
sum  of  powers  of  the  correlation  coefficient  and  it  is  assumed  that 
decorrelation  of  surface  heights  implies  statistical  independence. 

Using  these  joint  probability  density  functions,  Beckmann (197 3a, b) 
derives  physical  optics  and  geometrical  optics  approximations  for  the 
scattering  cross  sections.  Recently,  Brown  considered  rough  surfaces 
with  exponential  joint  height  probability  densities  (Brown,  1982) . 

In  his  work  Brown  concludes  that  the  appearance  of  a  delta  function  type 
specular  term  in  his  results  "suggests  that  decorrelation  does  indeed 
imply  statistical  independence  for  real  surfaces." 

In  Section  2  the  principal  elements  of  the  full  wave  approach  are 
summarized.  In  Sections  3  and  4,  the  expressions  for  the  marginal  slope 
densities  and  the  shadow  functions  for  the  non-Gaussian  surfaces  are 
derived.  In  Section  5  Illustrative  examples  are  presented.  By  considering 
a  very  broad  family  of  non-Gaussian  rough  surfaces  that  include  in  the 
limits  the  exponential  and  the  Gaussian  surfaces,  it  is  possible  to 
derive  the  full  wave  solutions  for  the  backscatter  cross  sections  for 
more  realistic  models  of  rough  surfaces.  Since  it  is  assumed  in  this 
work  that  decorrelation  of  the  surface  heights  implies  statistical 
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independence,  no  delta  function  type  specular  terns  are  obtained  in 
these  results.  It  is  found  that  for  angles  of  Incidence  6*  >  25°,  the 
backscatter  cross  section  for  the  horizontally  polarized  waves  are 
more  sensitive  to  the  specific  form  of  the  joint  height  probability 
density  assumed.  The  reason  for.  this  is  given  in  Section  5.  For  near 
normal  incidence  both  the  backscatter  cross  sections  for  vertically  and 
horizontally  polarized  waves,  exhibit  the  same  dependence  on  the  special 
form  of  the  joint  height  probability  density  assumed.  The  depolarized 
backscatter  cross  sections  are  sensitive  to  the  special  form  of  the 
surface  statistics  only  near  normal  Incidence. 

It  is  interesting  to  note  that  for  the  mean  square  slopes  considered 
in  the  illustrative  examples,  the  shadow  functions  for  the  non-Gaussian 
surfaces  are  practically  indistinguishable  from  the  shadow  function 
for  the  Gaussian  surface,  nevertheless  the  backscatter  cross  sections 
are  sensitive  to  the  precise  form  of  the  surface  height  statistics 
assumed  (see  Section  5). 


2.  Formulation  of  the  Problem 

The  full  wave  solutions  for  the  scattering  cross  sections  are  summarized 
in  this  section  for  a  two  scale  model  of  random  rough  surfaces  h(x,z)  (See 
Fig.  1)  (Bahar  et  al  1982).  To  this  end  consider  the  Fourier  transform  of 
the  surface  height  autocorrelation  function  <h(x,z)h' (x' ,z ')>  which  is  the 
surface  height  spectral  density  function  (Rice  1951) 


W(vx,vz)  -  -j  exp(ivxxd+  iv2zd)<h(x,z)h'(x’,z,)>dxddzd 
IT 


(2.1) 


In  (2.1)  the  symbol  <  > denotes  statistical  average  and 

v  ■  k  (n^-n*)  -  va  +  va  +  va  ,  Ivl  ■  v 
o  xx  yy  zz1 


(2.2a) 


in  which  n*  and  n^  are  unit  vectors  in  the  directions  of  the  incident  and 


scattered  waves 


n1  -  sinB1"  cos^1  a  -  cose1  a  +  sinO*  sintfi*  a  ■  kVk 
o  x  O  y  o  z  o 

n^  ■  sin6^  cos<j>^  a  +  cos8^  a  +  sln8^  sin<^  a  ■  k^/k 
o  x  oy  o  z  o 


(2.2b) 


(2.2c) 


and  k  ■  w(e  u  )  is  the  free  space  wavenumber  characterized  by  the  permittivity 
o  o  o 

and  permeability  An  exp(iut)  time  dependence  is  assumed  in  this  work. 

The  distance  between  two  points  in  the  reference  (x,z)  plane  is 


r.  «  |r.  ■  Ix.a  +  z.a  j  ■  |(x-x')a  +  (z-z')a 

d  1  d  1  d  x  d  z '  1 '  x  z 


(2.3) 


For  convenience  the  rough  surface  height  h(x,z)  is  decomposed  into  two  parts 

h(x,z)  -  h.(x,z)  +  h  (x,z)  (2.4) 

X  s 

the  first  term  h^  consists  of  the  large  scale  spectral  components 

0  <  k2  H  v2  +  v2  <  kj  (2.5) 

where  kd  is  the  wavenumber  where  spectral  splitting  is  assumed  to  occur 
(Brown  1978,  1980).  In  this  work  the  wavenumber  kd  is  chosen  such  that 
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the  surface  meets  the  radii  of  curvature  criteria  (associated  with  the 

Kirchhoff  approximations  for  the  surface  fields).  It  is  also  assumed  here 

that  the  surface  h^  meets  the  condition  for  deep  phase  modulation  (Bahar 

et  al  1982).  The  second  term  in  (2.4),  h  ,  consists  of  the  small  scale 

s 

spectral  components 

k .  <  k  <  k  (2.6) 

a  -  c 

where  k£  is  the  spectral  cutoff  wavenumber  (Brown  1978) .  The  full  wave 
scattering  cross  section  for  the  composite  surface  (2.4)  which  accounts  for 
both  specular  scatter  and  Bragg  scatter  in  a  self-consistent  manner  is  given 
by  the  weighted  sum  of  two  cross  sections  (Bahar  et  al  1982) 


<oP<^>T“  <oPQ>0  +  <<y*^>R  (2.7) 

in  which  the  first  and  second  superscripts  (P  and  Q,  V  for  vertical,  H  for 
horizontal)  denote  the  polarization  of  the  scattered  and  incident  waves 
respectively.  The  first  term  in  (2.7)  is  the  scattering  cross  section 


associated  with  the  surface  h^.  Thus 


<oPQ>«  ■  !xs(v)|2  <0PQ> 


where 


PQ  47Tko 
<0  y>  * 

CO 

V 

y 


n*a 


P_(nf  ,ni|n)p(ii) 


(2.8) 

(2.9) 


n+n 


In  (2.8)  the  characteristic  function  for  the  surface  h  ,  Y  (v) ,  has  the  effect 

s 

of  decreasing  the  contribution  of  the  physical  optics  scattering  cross 
PQ  2 

section  <0^ >.  As  <h  >,  the  mean  square  of  the  small  scale  surface  height 

.  s,2 

(that  rides  on  the  large  scale  surface  h^) ,  approaches  zero,  |x  |  ♦l. 

It  is  assumed  in  this  work  that  the  small  scale  surface  height,  h8,  has  a 
Gaussian  distribution  and  in  this  case  the  coefficient  |x*|  in  (2.8)  is 
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(2.10) 


|x8(v)|2  -  (<exp  iv  hg>|2  -  exp(-4k2<hg>)  -  exp(-g) 

I  8  l  ^ 

Thus  Che  effective  weighting  function,  |x  [  »  decreases  monotonlcally  as 

22  _f  _i  _ 

0  -  4kQ<hg>  increases.  The  shadow  function  Pj(n  ,n  |n)  is  the  probability 
that  a  point  on  the  rough  surface  is  both  illuminated  by  the  source  and 
visible  at  the  observation  point  given  the  value  of  the  slopes  hx»hz  at 
that  point  (Smith  1967 ,  Sancer  1969) .  The  unit  vector  normal  to  the  rough 
surface  is 

n  -  V(y-hA(x,z)}/ jvfy-h^x.z))  |  =  7(y-h(x,z))/|v(  y-h(x,z))  | 

-  (-hax+  ay-  hzaz)/(hx  +  1  +  h*)**  (2.11a) 

where 

h  ■  3h/3x  ,  h  -  3h/3z  (2.11b) 

X  Z 

and 

n  *  v/v  (2.11c) 

8 

is  the  value  of  n  at  the  specular  points  of  the  rough  surface.  Furthermore, 
p(n)  *  p(h  ,h  )  is  the  joint  probability  density  function  for  the  slopes  of 
the  large  scale  surface  h^  .  In  this  work  the  joint  probability  density  func¬ 
tion  for  the  large  scale  surface  height  h^(x,z)  and  its  associated  slope 
distribution  function  are  assumed  to  be  either  Gaussian  or  non-Gaussian 
(Beckmann  1973a, b,  see  Section  3). 

The  coefficient  in  (2.9)  depends  on  the  polarization  of  the  scattered 
and  incident  waves,  the  unit  vectors  n*,n^  and  n  and  the  permittivity  e  and 
permeability  y  of  the  media  above  and  below  the  rough  surface  h(x,z) 

(Bahar  1981a, b). 

The  second  term  in  (2.7)  is  the  scattering  cross  section  associated 

with  the  small  scale  surface  h  (that  rides  on  the  large  scale  surface). 

s 


Thus  (Bahar  et  al,  1982) 


|DPV(v-)|2  fv;l2m  «B(Vj,vf) 


<0PQ>  .  4k2  J  4 


au 

•P2(iif ,ni|ii)p(n)dhxdhz  -  £  <aP(^>Rm 

m«l 

In  (2.13)  V  is  the  two  dimensional  Fourier  transform 

IQ 

j"2L~z"  “  ^2  f  <hi hi>m  exp(1Vd+  iVd)dxdd?d 

■  |  Vi(vvz:)wi{vvvvf)dvx^ 

"iWm-l(vx‘V  ©  W1{VV 


(2.13) 


(2.13) 


(2.14) 


where  the  symbol  (x)  denotes  the  two  dimensional  convolution  of  with 

W^,  the  surface  height  spectral  density  function  for  the  small  scale  surface 


h  .  Thus 
s 


WV 


V(v-,Vr.)  for  k  >  k , 
x  z  d 


for  k  <  k^ 


(2.15) 


and  v-  and  v?  are  the  orthogonal  components  of  the  vector  v  (2.2)  tangent 

to  the  surface  h.(x,z)  while  v-  is  the  component  of  v  normal  to  the  surface  h.. 

)L  y  * 

Thus  v  in  the  local  coordinate  system  associated  with  the  rough  surface  h.  can 


be  expressed  as  (see  Fig.  2) 


where 


v  *  v  n,  +  v  n_  +  v  n» 
x  1  y  2  z  3 


(2.16a) 


nx  -  (n  x  az)/|n  xa'J.yn,  n3»  ^x  n 


(2.16b) 


VX  ”  v*nl  *  vy  "  v’n2  *  v*  "  v*n3 


(2.16c) 


f 


Using  the  full  wave  approach  it  is  shown  that  a  suitable  value  for  k.  (the 

a 

wavenumber  where  spectral  splitting  is  assumed  to  occur)  is  obtained  by 
2  2 

setting  0  ■  4k  <h  >-1.0  since  for  this  value  of  0  the  assumed  condition 
o  s 

for  deep  phase  modulation  is  satisfied.  For  0  >  1.0  the  total  scattering 

cross  section  (2.7)  is  insensitive  to  variations  in  the  value  of  k,  (Bahar 

a 

et  al  1982).  However,  since  k^  decreases  as  0  increases,  it  is  necessary 

to  evaluate  more  terms  in  (2.13)  for  larger  values  of  0. 

3.  Non-Gaussian  Joint  Probability  Density  Functions  for  the 
Large  Scale  Surface  Heights  and  the  Associated  Slope 
Distribution  Function 

The  principal  purpose  of  this  work  is  to  apply  the  full  wave  appi 
to  scattering  from  rough  surfaces  to  surfaces  with  non-Gaussian  surfa< 
height  distribution  functions  and  to  compare  the  results  with  these  for 
Gaussian  surfaces.  In  his  recent  work,  Brown  (1982)  investigates  scattering 
by  surfaces  for  which  decorrelation  does  not  imply  statistical  independence. 
In  this  work  Brown  concludes  that  the  appearance  of  a  specular  term  coupled 
with  the  lack  of  any  experimental  scattering  data  supporting  this  result 
"suggests  that  decorrelation  does  indeed  imply  statistical  independence  for 
real  surfaces . " 

Thus,  even  though  from  a  mathematical  point  of  view,  statistical 
independence  is  not,  in  general,  implied  by  a  lack  of  correlation  (Papoulis 
1965),  it  is  assumed  in  this  work  that  the  joint  surface  height  distribution 
function  becomes  statistically  independent  as  the  surface  height  decorrelates. 
Beckmann  (1973a, b)  derives  a  method  for  finding  the  joint  distribution 
function  for  the  surface  height,  f(h,h';rd),  from  a  given  marginal  p(h) 
that  is  not  necessarily  normal  and  s  given  correlation  coefficient  R(rd). 
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rr  *  Rn(r.) 

<hh’>  -  hh'f(h,h';rd)dhdh'  -  £  - - 

*  *  n*0 

-  R(rd)(<h2>  -  <h>2)  +  <h>2  (3.8) 

in  agreement:  with  the  condition  f(h,h';rd>  must  fulfill  if  it  is  to  be  a 
two  dimensional  density  of  a  random  process  with  given  correlation  coefficient 
R(rd).  In  addition  (3.1)  clearly  satisfies 


Aim  f(h,h';r.)  -*■  p(h)p(h')  (3.9) 

r.-Hs  d 

Q 

Finally  express  the  function  g(h)  as  follows  in  an  infinite  series  of 
orthogonal  polynomials 

00 

g(h)  -  l  On  Qn(h)  (3.10) 

n-0 

where  on  using  (3.4)  the  constant  is 


g(h')p(h')Qn(h,)dh’/qn 


On  interchanging  integration  on  summation  (3.10)  can  be  written  as 


g(h) 


g(h')[p(h’)  l  ~  Qn(h)Qn(h,)]dhI 
n-0  qn 


(3.11) 


(3.12) 


thus  the  quantity  in  square  brackets  is  the  Dirac  delta  function  6(h-h'). 
Thus  from  (3.1)  it  follows  that 


Aim  f(h,h';r.)  *  p(h)6(h-h')  ■  p(h')6(h-h')  *(3.13) 

r  -*-o  d 

d 

On  integrating  (3.13)  with  respect  to  h  and  h'  it  is  readily  shown  to  satisfy 

(3,3).  (Contrary  to  the  statements  by  Beckmann  (1973a, b),  the  Dirac  delta 

function  does  not  remain  unchanged  when  multiplied  by  p(h)  and  Aim  f(h,h';r.) 

-  r^o  a 

+  6(h-h')j.  In  Beckmann's  work  (1973a, b),  the  orthogonal  polynomials  (^(h) 
can  be  the  Hermite  polynomials,  the  generalized  Laguerre  polynomials  or  the 
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Jacobi  polynomials.  However,  in  this  work  it  is  assumed  for  convenience 
that  are  the  generalized  Laguerre  polynomials  (Abramowitz  and  Stegun 
(1964)) 

Q^h)  -  L^k_1)(h)  (3.14) 

Since  the  weighting  functions  associated  with  the  Laguerre  polynomials  are 
proportional  to  the  gamma  functions,  the  marginals  for  the  large  scale 
surface  height  probability  density  functions  are  the  gamma  distributions 

k  k-1  -vh 

p(h)  -  (kIl)-> -  ,  h  >  0  (3.15) 

2  2 

for  which  the  mean  is  <h>  *  k/v,  the  variance  is o  ■  k/v  and  the  normalization 
2 

constant  (3.4)  is 

q^  -  (k  +  n-l).7(k-l).'n.'  *(3.16) 

For  k  *  1  (3.15)  reduces  to  the  exponential  distribution 

p(h)  »  ve  Vh  ,  (k*l)  (3.17) 

and  for  the  limiting  case  k  ®»  (3.15)  reduces  to  the  Gaussian  distribution 
by  virtue  of  the  Central  Limit  Theorem.  (For  most  practical  purposes 
k  >  25  is  sufficiently  close  to  the  limiting  case). 


p(h)  -  -  exp  -[ - =- — ],(k*«)  (3. 

a/B  2a 

2 

where  a  is  the  variance  and  <h>  is  the  mean  height. 

Thus  using  the  gamma  distribution  (3.15),  it  is  sufficient  to  vary  one 
parameter  k  to  examine  the  effects  of  changing  the  marginal  height  distribu¬ 
tion  from  exponential  to  normal  in  gradual  steps. 

Following  the  procedures  outlined  above,  the  joint  distribution  for 


(3.18) 


the  surface  heights  whose  marginals  are  p(h)  (3.15)  (with  v  ■  1  and  variance 
k  can  be  shown  to  be  given  by  (Beckmann  1973b)) 


1 


,  (hh')k-1  e_(h+h  J  r  R  CrdJn’  k-1,  .  k-1.. 

f(h.h’;r.)  -  '  -  I  -nrS3VrLn  <h>Ln  <h  > 


[(k-1)::  »:o  <k+”-1):  " 

k-1 

in  which  L  (h)  are  Laguerre  orthogonal  polynomials  (Abramovitz  and  Stegun 
n 

1964).  The  corresponding  (one  dimensional)  slope  distribution  function 


(3.19) 


p(hx)  can  be  obtained  from  the  inverse  Fourier  transform 


p(hx> 


r j  e  »  R  <rj)  r 

~7~  1 


p(h)Lk-1(h)exp(iwh) |2  exp(-irdhxw)dw 


where 


(3.20) 


oo  n 

I  h 

n-o  qn 


p<h)Lk-1(h)«pa»k)dhj-j 

1  n-0  n:  (1+w  )  (k-1) 


[1  +  (1-R)v2]k 


*(3.21) 


Thus 


p(h  )  -  Urn 

x  r.-*-o 
a 


exp(-iwhx)dw 


~~  2.  2.k 
(1  +  w  /w  ) 
o 


(3.22a) 


.  th 


where  the  k  order  poles  are  at  points 


w  -  +  1/SI-R  •  t  i  w  (3.22b) 

±  -  o 


Using  standard  procedures  for  integration  in  the  complex  w  plane,  the  path 
along  the  real  axis  is  closed  by  an  infinite  semicircle  in  the  lower  or 
upper  half  plane  depending  on  whether  z  >  0  or  z  <  0.  Thus  the  pole 


enclosed  by  the  contour  is  at 

w  »**i  wq,  for  z  >  0 

w  -  *(3.23) 

w  ■  i  w  ,  for  z  <  0 
+  o 

and  (3.22)  reduces  to 
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*(3.24) 


Ye"Ylhxl  k-1  (2k-j-2)!(2Y|hx|)J 

p(h  )  ”  2k-l  £  1  * (k-1- 1) ! 

X  2Zk  1 (k-1) I  j-0  3  '  3  ' 


where 


Y  -  Aim  r.  w  -  Aim  r.//l-R 

r<j-hj  d  °  “ 


(3.25) 


For  a  stationary  random  process,  the  correlation  function  is  an  even  function 


of  r^  (Beckmann  1968)  and  for  small  r^ 

_2 


*<V  ■  1  -  ^2  rd  +  °<^> 


(3.26) 


where  is  the  mean  square  slope 


a  -  -a*  R' '  (o) 
x 

Therefore  (3.25)  reduces  to 

ST  a  /2k 


(3.27) 


(v-1) 


*(3.28) 


x  x 

and  the  slope  distribution  function  corresponding  to  the  two  dimensional 
gamma  (surface  height)  distribution  (3.19)  is 
2/2k 


P(hx) 


exp(-/2klh  1  /o  )  k-1  ,  . 


*(3.29) 


The  above  result  can  be  shown  to  hold  for  arbitrary  values  of  v  and  not  just 
v-1  (Beckmann  1973a).  It  is  assumed  in  this  work  that  the  probability 


density  for  the  slopes  hx  and  h^  are  independent  (as  for  the  Gaussian  case 


Barrick  1970,  Brown  1978,  Bahar  1981a, b) ,  thus  the  two  dimensional  slope 
distribution  for  the  large  scale  surface  is 

p(n)  -  p(h  ,h  )  -  p(h  )p(h  )  (3.30) 

X  Z  x  z 


in  which  p(h  )  and  p(h  )  are  the  same  functions  as  in  (3.29).  The  corresponding 
X  z 


slope  distribution  function  is 


P(bx) 


«*p(”bV2o  ) 

Si T  a 


(3.31) 
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4.  Th«  Explicit  Expression  for  Cho  Shadow  Functions  Associated 
with  the  Mon-Csusslan  Surface  Height  Oietrlhition  Functions 


In  this  section  explici'  -ess ions  are  derived  for  the  shadow  functions 
related  to  the  non-Gaussian  surface  height  distribution  functions  considered 
in  Section  3.  Brown  (1980)  recently  developed  a  very  convenient  procedure 
to  evaluate  the  shadow  functions  for  surfaces  with  non-Gaussian  height  distrib¬ 
ution.  Following  Smith’s  work  (1967),  Brown  (1980)  assumes  that  the  surface 
heights  h  and  slopes  hx  at  points  separated  by  the  distance  are  uncorrelated. 
In  addition  Brown  assumes  that  decorrelation  implies  statistical  Independence. 
Thus  keeping  the  notation  used  in  this  paper,  the  shadow  function  for  back- 
scatter  (n1-  -n1),  ?2 (-n^n1 |n) ,  is  the  probability  that  a  point  on  the 
rough  surface  is  illuminated  by  the  source  (and  visible  at  the  observation 
point)  given  the  value  of  the  slope  at  that  point. 

P2(-n±.ni|S)  -  S(-n1*n)/[l  +  rVy1]  (4.1) 

where  S(-n1*n)  is  the  unit  step  function.  The  argument  -ni*n  vanishes  when 
the  incident  (or  backscattered)  wave  normal  is  tangent  to  the  surface. 

Thus,  when  the  plane  of  incidence  is  the  x,y  plane  (see  Fig.  1),  the  argument 
(-n^n)  vanishes  for 

y1  =  cotS1  ■  3h/3x  =  hx  (4.2) 

in  which  81  is  the  angle  of  incidence  (2.26)  with  respect  to  the  reference 
surface  (the  x,z  plane).  Furthermore,  in  (4.1) 


r1  s  r^1)  -  |  (hx-yi)P(hx)dhx 
y1 


(4.3) 


in  which  p(hx)  is  the  one  dimensional  slope  distribution  function.  On 
substituting  the  non-Gaussian  slope  distribution  function  (3.29)  (associated 
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with  the  gamma  surface  height  distribution)  it  can  be  shown  that 


iw  i  „  exp(-^  kr^  (2k-1-2).'  2J  jrX  r(l^  ^qx)J 

''u  "  2k-l  i  (k-1-1)  ’  “  (1+1-r) 1 

2^k  (k-l)!  j-0  J  L)’  r-1  0+1  1,1 


T(y  )/y 


The  shadow  function  for  the  bistatic  case,  n%  -n*,  is  (Brown  1980) 


P2(nf  .n^n) 


SO-n^n)  Qi  ^  0f  xf  xi _ 

- i — 7~  •  e„  >  e«  »  9  -9  *  11 

i+ri/yi  0  “  ° 

S(nf*n)  „f  .  „i  ^f  ,i 

— r~r  •  e0  -  eo  *  ♦  -  *  • n 


1+r/y 

S(n^*n)S(-n^*n) 


i+(ri/yi)+(rf/yf) 


elsewhere 


(4.5) 


In  (4.5) 


Tf  =  T(yf)  and  yf  -  cot8f  (4.6) 

o 

where  0^  is  the  scatter  angle  (2.2c)  with  respect  to  the  reference  plane. 
The  above  expressions  (4.5)  are  also  in  agreement  with  Sancer's  (1969) 
results  for  Gaussian  surface  height  distributions.  (In  Sancer's  published 
results  (1969),  the  inequality  symbols  should  be  reversed,  Brown  (1980)) . 


5.  Illustrative  Examples 

In  this  section  the  like  and  depolarized  backscatter  cross  sections 
are  evaluated  for  several  composite  models  of  perfectly  conducting  rough 
surfaces  with  different  surface  height  distribution  functions.  The  surface 
h^  consisting  of  the  large  scale  spectral  components,  0  <  k  <  kj,  is  assumed 
to  be  characterized  by  the  ganoa  surface  height  distribution  function  (3.15) 
and  the  associated  non-Gaussian  slope  distribution  function  (3.29).  For 
the  illustrative  examples  presented  here  the  parameter  k  is  set  equal  to 
1,2,5  and  25  and  the  limiting  case  (k  °°)  is  represented  by  the  Gaussian 
distribution  (3.18).  Two  different  total  mean  square  slopes  are  considered 
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2  2 

(a)  a^t  -  0.0171  and  (b)  -  0.0564.  The  surface  hg  consisting  of  the 

small  scale  spectral  components,  k  >  k. ,  is  characterized  by  Gaussian 

—  a 

surface  height  distribution  functions.  The  wavenumber  kd  where  spectral 

2  2 

splitting  is  assumed  to  occur  is  determined  by  setting  8  ■  4kQ<hg>  ■  1 
(Bahar  et  al  1982)  and  the  specific  form  assumed  for  the  surface  height 
spectral  density  function  (2.1)  is  (Brown  1978) 


W(k)  -  ~  S(k) 


2  2  4 

(k%  o* 


,  k  <  k 
-  c 


I  0  ,  k  >  k  (5.1) 

v  c 

In  (5.1)  w(k)  is  the  spectral  density  function  employed  by  Rice  (1951) 
and 


B  *  .0046  (case  a)  ;  B  ■  0.0133  (case  b) 

k2  -  v2  +  v-  (cm)”^  ;  k  -  12  (cm)  1 

X  z  c 

K  -  (335.2  V4)"1*  (cm)-1  ;  V  -  4.3  (m/s)  (5.2) 

The  wavelength  for  the  electromagnetic  wave  is 

\q  -  2  (cm)  ,  (kQ  -  3.1416  (cm)"1)  (5.3) 

The  mean  square  height  for  the  small  scale  surface  hg  is 
2it  k_ 

(5.4) 
d 


<h2> 

s 


W(k) 

4 


kdkd4>  - 


and  the  total  mean  square  slope  for  the  large  scale  surface  h^  is  (Brown 
1978) 


IT  k. 


2  2  ,  2 

o„  ■  a  +  a  ■ 

It  x  z 


k3  dkd<p  = 


- 

~~ 

2  2 
kZ  +  K 

d 

L  K  J 

(5.5) 


o  o 


2  2 

in  which  0  ■  O  ■ 

X  z 
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In  Fig.  3,  a  times  the  one  dimensional  ganma  surface  height  distribu¬ 
tions  (3.15)  are  plotted  as  functions  of  h/o  for  k  -  1,2,5,25  together  with 
the  Gaussian  surface  height  distribution  (with  <h>/0  ■  &  m  5) .  In  Fig.  4, 
the  corresponding  one  dimensional  slope  distributions  times  0^x  are 

plotted  as  functions  of  h  /o^.  In  Fi8*  5  the  corresponding  shadow 

functions  (4.5)  are  plotted  as  functions  of  6*  for  backscatter  (n^»  -nS 
2  2 

(a)  at  “  0.0171,  (b)  o^t  ■  0.0564.  As  expected  the  shadow  function  is 

smaller  for  the  larger  mean  slopes.  The  shadow  functions  are  insensitive 

to  the  particular  form  of  the  slope  density  functions  considered. 

In  Figs.  6a,b,c  and  7a,b,c,  the  corresponding  like  polarized  and 

depolarized  total  backscatter  cross  sections  (2.7)  are  plotted  as  functions 

of  6i  (a)  <oVV>  (b)  <0HH>  (c)  <oVH>  *  <o^V>  in  Figs.  6a,b,c  o?  «  0.0171 

o  xt 

2 

and  in  Figs.  7a,b,c  «  0.0564.  The  like  backscatter  cross  section 

near  normal  incidence  (6  -  0)  is  largest  for  the  case  k  *  1.  This  is 

o 

expected  since  the  case  k  *  1  corresponds  to  a  surface  with  a  slope 
distribution  which  has  the  largest  zero  slope  probability.  On  the  other 
hand  consistent  with  the  above  results,  the  depolarized  backscatter  cross 
sections  near  normal  incidence  are  smallest  for  the  case  k  ■  1.  The  back¬ 
scatter  cross  section  for  k  «  25  is  practically  Indistinguishable  from 
the  backscatter  for  Gaussian  random  surfaces.  Furthermore,  as  one  may 
expect,  the  spread  in  the  values  for  the  backscatter  cross  sections  as  k 
varies  from  1  through  00  increases  as  the  mean  square  slope  0£t  increases. 

At  large  angles  of  incidence  0*,  the  spread  in  the  values  of  the  like 
backscatter  cross  sections  is  more  pronounced  for  the  horizontally  polarized 
waves  than  for  the  vertically  polarized  waves.  The  reason  for  this  is  found 


by  examining  the  expression  £or  the  dominant  term  at  large  angles  0^,  i.e., 
<oPQ>ri  (first  order  Bragg).  In  this  term  |v^0^|^  is  proportional  to 
In1  while  Iv-D^l^  is  proportional  to  |2-(n*  *n)^  The  latter  term  Is 

both  larger  and  less  sensitive  to  slope  variations  n(hx,hz).  In  Figs.  8a,b,c,d 

VV  W 

and  9a,b,c,d,<0  >Tthetotal  backscatter  cross  sections, <o  >q  the  large  scale 

W 

backscatter  cross  sectioned  t he  leading  term  of  the  small  scale  cross 

section  (which  corresponds  to  first  order  Bragg  scatter)  and  the  remainder 

term<a^^+<CT^>^are  plotted  as  functions  of  0^,  case  (a)  k  *  1,  case  (b) 

k  *  2,  case  (c)  k  -  5,  case  (d)  Gaussian  surface  heights.  The  results  for 

k  *  25  are  not  shown  since  they  are  indistinguishable  from  the  results  for 

2 

the  Gaussian  surface).  In  Figs.  8a,b,c,d,  a ^  «  0.0171  and  in  Figs. 

2  HH  HV  VH 

9a,b,c,d  o^t  -  0.0564.  The  corresponding  results  for  o  and  o  ■  o 

are  presented  in  Figs.  10a,b,c,d  and  lla,b,c,d,  12a,b,c,d  and  13a,b,c,d. 

From  the  above  plots  of  the  backscatter  cross  section  it  is  seen  that  not 

only  is  the  total  backscatter  cross  sect ion <  >T  sensitive  to  changes  in 

values  of  k,  but  the  individual  terms  in  the  weighted  sum  (2.7)  are  also 

sensitive  to  changes  in  k.  Furthermore,  except  for  near  grazing  angles 

PQ 

the  backscatter  cross  section  for  the  small  scale  surface  <o  cannot 
be  approximated  by  the  leading  term  in  the  series  (2.13)  which  corresponds 
to  first  order  Bragg  scatter.  Ic  has  been  shown  using  the  full  wave 
approach,  Bahar  et  al.  (1982),  that  the  total  backscatter  cross  section 
<0P%  is  Independent  of  the  specific  choice  of  kj  for  8  >  1  provided  none 
of  the  significant  terms  in  (2.13)  are  neglected.  The  terms  <tp  , for  m  >  4, 
can  be  neglected  for  all  0q  when  k^  is  chosen  such  that  8-1. 

On  the  basis  of  the  above  results,  several  schemes  can  be  devised  to 
distinguish  between  the  backscatter  cross  sections  for  different  values  of  k. 


The  behavior  of  «^P>T  (P-V.H)  near  normal  incidence  is  of  particular 


interest  in  this  respect.  Thus  one  notices  that  |d<aPVd0*|  i  is  largest 


for  k  *  1  and  decreases  gradually  as  k  increases. 


UU 

Since  <0  is 


more  sensitive  to  variations  in  k  for  6  >  25  ,  the  backscatter  cross 

o 


sections  for  horizontally  polarized  waves  are  more  indicative  of  the 
particular  form  of  the  rough  surfaces  height  distribution  than  the  backscatter 
cross  sections  for  the  vertically  polarized  waves.  The  depolarized 

backscatter  cross  sections  <0  >(Pi*Q)  are  rather  insensitive  to  the  particular 

form  of  the  rough  surface  height  distribution  except  near  normal  incidence 


where 


|d<aP<fyd6c 


IV0 


is  largest  for  k  -  1. 


6.  Concluding  Remarks 

Full  wave  expressions  for  the  like  polarized  and  the  depolarized 
scattering  cross  sections  are  derived  for  a  broad  family  of  non-Gaussian 
rough  surfaces.  It  is  assumed  that  when  the  surface  heights  decorrelate 
they  become  statistically  independent.  For  the  cases  considered  in  the 
illustrative  examples,  it  is  shown  that  while  the  shadow  functions  are  not 
very  sensitive  to  the  parameter  k(*l,2,3...)  of  the  surface  height 
probability  density  assumed,  the  backscatter  cross  sections  are  sensitive 
to  the  parameter  k.  The  like  polarized  back- 

vv  hh 

scatter  cross  sections  <o  >T  and  <o  >T  have  the  same  dependence  on  k 

for  small  angles  of  incidence  0*,  however  for  large  angles  0*,  the  cross 

HH  VV 

section  <o  is  much  more  sensitive  to  k  than  <a  >^.  An  examination  of 

the  expression  for  the  dominant  term  corresponding  to  first  order  Bragg 
PP 

scatter  <o  >R1  provides  an  explanation  to  this  phenomenon.  The  depolarized 

PQ 

backscatter  cross  sections  <0  >T  (P/Q)  are  sensitive  to  k  only  near  normal 
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Incidence.  For  8^  -  0,  Che  like  polarized  cross  section  <c  >T  is  largest 

for  the  smallest  value  of  k  (k  •  1).  This  is  because  the  slope  probability 

density  for  zero  slope  is  largest  when  k  -  1.  Thus  consistent  with  energy 

conservation  as  the  like  polarized  cross  section  increases  with  decreasing 

k  values,  the  depolarized  cross  section  decreases  with  decreasing  k  values. 

Both  the  total  backscatter  cross  sections  and  the  components  of  the  weighted 

sum  of  cross  sections  are  presented  in  Section  5.  The  terms  corresponding 
PQ 

to  <0  m  >  2  cannot  be  neglected  in  this  work  and  a  perturbed  physical 

optics  approach  is  not  suitable  for  determining  the  total  backscatter 
cross  section  (Bahar  et  al.  1982). 
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Figure  Captions 


Figure  1 
Figure  2 

Figure  3 

Figure  4 

Figure  5 

Figure  6 

Figure  7 

Figure  8. 

Figure  9. 


Plane  of  incidence,  scattering  plane,  and  reference  (x.z)  plane. 

Local  plane  of  incidence  and  acatter  and  local  coordinate  system 
with  unit  vectors  n^.n^.n^. 

The  one  dimensional  gamma  surface  height  distributions  times  o 
as  functions  of  h/o  for  k  -  1,2,5,25  together  with  the  Gaussian 
surface  height  distribution  with  <h>/o  «  5. 

(O)  k  -  1,  (A)  k  *  2,  (+)  k  -  5,  (X)  k  -  25,  (O)  Gaussian. 

The  one  dimensional  slope  distributions  times  *  /l  as 
functions  of  (O)  k  »  1,  (A)  k  »  2,  (+)  k  •  5,  (X)  k  *  25, 

(O0  Gaussian. 


The  shadow  functions  of  functions  of  6q  for  backscatter. 

a\t  -  0.0171  and  -  0.0564  Q  k  -  1,  (A)  k  -  2. 

(+)  k  -  5,  (X)  k  -  25,  (O)  Gaussian. 

The  total  backscatter  cross  sections  (2.7)  as  functions  of  61 


for  ct  -  0.0171 
XtW 

(a)  <a  > 


(b)  <oHH> 


(c)  «,"*>  -<**>. 


(O  k  -  1,  (A)  k  -  2,  (+)  k  -  5,  (X)  k  -  25,  (0)  Gaussian. 

The  total  backscatter  cross  sections  (2.7)  as  functions  of  0*  for 


(Jlt-  0.0564 

vv 

(a)  <o  > 


(b)  <aHH> 


(c)  <aVH>  ■  <oHV> 

(O)  k  -  1,  (A)  k  -  2,  (+)  k  -  5,  (X)  k  -  25,  (<>)  Gaussian. 

Backscatter  cross  sections  (-^cr^^,  VO<^>0,  (D)<o^V^(A)<oVV^2+<oVV 


i  2 

as  functions  of  6  for  o„.  ■  0.0171 
o  it 


(a)  k  -  1, 


(b)  k  -  2, 


(c)  k  -  5, 


(d)  Gaussian 


W  W  VV  w  VV 

Backscatter  cross  sections  («)<a  >^,(X)<a  >,  (0)<a  ^,(A )<a  >Rj 

i  2 

as  functions  of  6  for  O.^  ■  0.0564 
o  It 


(a)  k  -  1, 


(b)  k  -  2, 


(c)  k  -  5, 


(d)  Gaussian. 


Figure  10.  Backscatter  cross  sections  (-)<c?\,  (X)<<J®V  .(EDccr®^  , 

”  o  R1 

^as  functions  of  8*  for  0^t  ■  0.0171 

(a)  k  -  1,  (b)  k  ■  2,  (c)  k  ■  5,  (d)  Gaussian 

Figure  11.  Backscatter  cross  sections  (Xjcc®^,  (P)<oHH> 

To  R 1 

(A)<a®^<a®^s  functions  of  8*  for  a^t  *  0.0564 
(a)  k  ■  1,  (b)  k-2,  (c)  k  »  5,  (d)  Gaussian. 

Figure  12.  Backscatter  cross  sections  (~)<0  \(*Xfl\(axa\1( 

i  O  R1 

(A)<  oVI^< °V1^as  functions  of  8*  for  oj^  *  0.0171 

(a)  k  *  1,  (b)  k  “  2,  (c)  k  *  5,  (d)  Gaussian. 

Figure  13.  Backscatter  cross  sections  (-)<  oVI^,  (X^o^, 

(A)<  0VT^2+<0VI^as  functions  of  9*  for  o^.  *  0.0564. 

(a)  k  *  1,  (b)  k  ■  2,  (c)  k  *  5,  (d)  Gaussian. 
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